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Abstract 

Crum's theorem in one-dimensional quantum mechanics asserts the existence of 
an associated Hamihonian system for any given Hamihonian with the complete set of 
eigenvalues and eigenfunctions. The associated system is iso-spectral to the original one 
except for the lowest energy state, which is deleted. A modification due to Krein-Adler 
provides algebraic construction of a new complete Hamiltonian system by deleting a 
finite number of energy levels. Here we present a discrete version of the modification 
based on the Crum's theorem for the 'discrete' quantum mechanics developed by two 
of the present authors. 

1 Introduction 

Crum's seminal paper of 1955 [1] has played an essential role in elucidating the structure of 
one-dimensional quantum mechanical systems in general and exactly solvable ones, in partic- 
ular. Throughout this paper, we mean 'exact solvability' in the Schrodinger picture, namely 
a quantum system is exactly solved when the complete set of eigenvalues and eigenfunctions 
are known. Many exactly solvable quantum mechanical Hamiltonians were constructed and 
investigated by combining shape invariance [2] and Crum's theorem [H |3] , or the factorisa- 
tion method [1] or the method of the so-called supersymmetric quantum mechanics [5]. It is 
interesting to note that most of these shape invariant systems are also solvable in the Heisen- 
berg picture [6]. Exactly solvable quantum mechanical systems of one and many degrees of 



freedom are not only important in their own right but also have fundamental applications 
in various disciplines of physics/mathematics, e.g. the Fokker-Planck equations [7] and their 
discretised version, birth and death processes [8], to name a few. 

Shape invariance is a sufficient condition for exactly solvable quantum mechanical sys- 
tems. The number of shape invariant systems, however, was quite limited; only about a 
dozen until the recent discovery [9l [T0|, [11] of the several types of infinitely many shape 
invariant Hamiltonians [121 [131 [HI [El [H] which led to the infinitely many exceptional La- 
guerre, Jacobi, Wilson and Askey- Wilson polynomials. Many methods were proposed to 
derive exactly solvable but non-shape invariant quantum mechanical systems from known 
shape invariant ones [IZl [HI [191 [201 [21] . (We apologise to those whose work we have missed.) 
Among them Krein-Adler's modification [18] of Crum's theorem is the most comprehensive 
way to generate infinitely many variants of exactly solvable Hamiltonians and their eigen- 
functions, starting from an exactly solvable one. The derived system is iso-spectral with 
the original one except that a finite number of energy levels are deleted. If the original 
system has polynomial eigenfunctions, as is usually the case, the derived systems have also 
poljTiomial eigenfunctions. By construction, these polynomials constitute a complete set of 
orthogonal functions. But they do not qualify to be called exceptional orthogonal polynomials 
[T2l [T3l [TH im [16] since some members of certain degrees are missing due to the deletion. 

The discrete quantum mechanics is a deformation of the ordinary quantum mechanics 
in the sense that the Schrodinger equation is a second order difference equation instead of 
differential. In the formulation of Odake and Sasaki [221 l23l [2^ [25] , the algebraic and ana- 
lytical structure of quantum mechanics as well as shape invariance and exact solvability are 
retained in the discrete version. The eigenfunctions of the exactly solvable one- dimensional 
discrete quantum mechanics are the Askey-scheme of hypergeometric orthogonal polyno- 
mials and their g- versions [2S1 [23 [2S1 [2S], e.g. the continuous Hahn, the Wilson and the 
Askey- Wilson polynomials. These examples are all shape invariant and they are also solvable 
in the Heisenberg picture [221 [S]- The dynamical symmetry algebra of these algebras are the 
Askey- Wilson algebras [SUl IM] and degenerations, which contain the g-oscillator algebra [32] • 
The discrete version of Crum's theorem was also established recently [331 [31] • 

In this paper we present the discrete quantum mechanics version of Adler's modification 
[18] of Crum's theorem. It allows to generate an infinite variety of exactly solvable discrete 
quantum Hamiltonian systems. The insight obtained from Crum's theorems and their mod- 
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ification, in the ordinary and the discrete quantum mechanics, is essential for the recent 
derivation of the infinite numbers of shape invariant systems and the new exceptional or- 
thogonal polynomials [121 [131 US]. We will discuss the main results, the specialisation to the 
cases of polynomial eigenf unctions and simplest example for various exactly solvable cases; 
first for the ordinary quantum mechanics and then for the discrete versions. The reason is 
two-fold; firstly to introduce appropriate notion and notation in the familiar cases of the 
ordinary quantum mechanics. Secondly we choose to reveal the underlying logical processes 
which are not easy to fathom in Adler's paper [18j or in Crum's original article [1]. As seen 
in the subsequent sections, the logical structures of the associated Hamiltonian systems and 
their modification by deletion of energy levels are shared by the ordinary and the discrete 
quantum mechanics. 

This paper is organised as follows. In section two, Adler's modification of Crum's theorem 
is recapitulated in appropriate notation for our purposes. The specialisation to the cases of 
polynomial eigenfunctions is discussed in some detail. Section three provides the discrete 
quantum mechanics version of the modification of Crum's theorem. Again the specialisation 
to the cases of polynomial eigenfunctions is mentioned. Appendix gives the simplest examples 
of the modified Hamiltonian systems obtained by deleting the lowest lying i excited states 
for various exactly solvable Hamiltonians. Appendix A provides three examples from the 
ordinary quantum mechanics, the harmonic oscillator, the radial oscillator, the Darboux- 
Poschl- Teller potential. Appendix B is for the four examples from the discrete quantum 
mechanics, the Hamiltonians of the Meixner-Pollaczek, the continuous Halm, the Wilson 
and the Askey- Wilson polynomials |221 125], which are known to reduce to the Hermite, the 
Laguerre and the Jacobi polynomials in certain limits, respectively. 

2 Ordinary Quantum Mechanics 

2.1 Adler's modification of Crum's theorem 

Let us start with a generic one- dimensional quantum mechanical (QM) system having dis- 
crete semi-infinite energy levels only: 



= So < £i < £2 < ■ ■ ■ ■ 
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Here we have chosen the constant part of the Hamiltonian so that the groundstate energy is 
zero. Then the Hamiltonian is positive semi- definite and can be factorised, 



/dW{x)Y d^W{x) . d 



« = p^ + f/M = p^+^_^j p = (2.2) 

= ^^def d dWjx) ^ d dW{x) 

' dx dx ^ dx dx 

Here a real and smooth function yV{x) G C°° is called a prepotential and it parametrises the 
groundstate wavefunction (t>Q{x), which has no node and can be chosen real and positive: 

Mx) = e"^^"^- (2.4) 

It is trivial to verify 

AM^) = ^ nMx) = 0. (2.5) 

In one dimension all the energy levels are non-degenerate. By construction all the eigenfunc- 
tions are square-integrable and orthogonal with each other and form a complete basis of the 
Hilbert space: 



%4>n{x) = Sn(t)n{x), U G Z+, (2.6) 

(t>n{x)*4>m{x)dx^hn5nm, < /i„ < oo, n,meZ+, (2.7) 

'XI 

where Z+ is the set of non-negative integers {0, 1,2, . . .}. It is well-known that the ra-th 
excited wavefunction 0„(x) has n nodes in the interior. For simplicity we choose all the 
eigenfunctions to be real. A few exactly solvable examples are given in Appendix. 

Let us choose a set of i distinct non- negative integers V =^ {di, ^2, ■ ■ ■ , di} C Z+, satis- 
fying the condition 

i 

Y[{m - dj) > 0, Vm e Z+. (2.8) 

This condition means that the set T> consists of several clusters, each containing an even 
number of contiguous integers 

dki,dki + 1, ■ ■ • ,dk2 rffcg, dfcg -I- 1, • • • , dfe4 ; dfcg, dfcg -I- 1, • • • , dfcg, ; • • • , (2.9) 

where rffej + l < c^feg, dk^ + l < dk^, ■ ■ • . If dk^ — for the lowest lying cluster, it could contain 
an even or odd number of contiguous integers. The set T> specifies the energy levels to be 
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deleted. This simply reflects the fact that no singularity arises when two neighbouring levels 
are deleted. In the ordinary QM, the zeros of the two neighbouring eigenfunctions (pj and 
4>j+i interlace with each other. This fact is essential for the non-singularity of the potential 
after deletion. See Adler's paper [18] for a proof. The situation is essentially the same in 
the discrete QM. However, in the dQM to be discussed in the subsequent section, due to the 
lack of general theorem, the interlacing of the zeros of the two neighbouring eigenfunctions 
(pj and 4>j+i must be verified for each specific Hamiltonian. Deleting an arbitrary number of 
contiguous energy levels starting from the groundstate [dk^ = 0) is achieved by the original 
Crum's theorem [1]. 

Next we will construct Hamiltonian systems corresponding to the successive deletions 
Hdi... (and Ad^..., , etc.) step by step, algebraically. It should be noted that some 
quantities in the intermediate steps could be singular. For given di the first Hamiltonian 7i 
can be expressed in two different ways: 



n = AU = Al^Ad, + Sd„ Ad,<pd, = 0, (2.10) 
dcf d dWd,{x) t dcf d dWd^x) . . s dcf , , . . „ . 

" d^' ^'d^ d^' ^^^^^ " log</>di(^), (2-11) 

, _ f dWjx) ^^ , d'Wjx) _ f dWdAx) \' , d'WdAx) , ^ 

^"l dx J ^ dx^ ~\ dx J ^ dx^ + ^ ' 

Unless di = 0, Wdiix) is singular due to the zeros of (pdiix). It is very important to note 
that in (|2TT]) is a 'formal adjoint' of Adi- We stick to this notation, since the algebraic 
structure of various expressions appearing in the deletion processes, from fl2.10p to fl2.34p . 
are best described by using the 'formal adjoint'. These define a new Hamiltonian system 

Tid, = Ad.Al^ + £d, = / + Ud,{x), (2.13) 

dcf 



Ud^ix) 

with the 'eigenfunctions' 



'Hdi4>din{x) = Sn(pdin{x), (pdi n{x) = Adi(pnix) , U E Z+\{di} . (2.15) 

Note that the energy level di is now deleted, (pdidi{x) = 0, from the set of 'eigenfunctions' 
{(pdinix)} of the new Hamiltonian Hdi- 

Suppose we have determined T-Ldi...ds and (pdi ...d^nix) with s deletions. They have the 
following properties 

nd,...d. =Ad,...dAi...d^+£d. =p' + Ud,...dM, (2.16) 



>Vrf,...d,(x)'=^'log0,,...,,(x), (2.17) 

. def d dWd,...d,{x) ^ def d dWd,...ds{x) , . 

- ^ j + Sd^, (2.19) 

0di...(i.n(a;) = ^di...<i>di...d,_in(2;) (n e Z+\{(^i, . . . ,4}), (2.20) 

'Hd^...ds(l>di...dsn{x) ^ Sn(t>di...dsn{x) (n G Z+\{(ii, . . . , dj) . (2.21) 

We have also 

4>d,...d,.,n{.x) = -A^^,^ _ ,^^(^) (n G Z+\K,...,4}). (2.22) 

^ds 

Next we will define a new Hamiltonian system with one more deletion of the level ds+i. We 
can show 

'Hdi...ds — '^l.i ...dsds+i'^di ...dsds+i ''<' ^ds+i, >^di ...ds ds+i<l>di ...da d^+i (2.23) 

. def d dW,i,...d,,L^Ax) def d d.W,i,...d,.d,+Ax) , 
Ad,...dada^,- , Ad,...dada^,- > {^-^V 

yVdi...dsds+i{x) = log0di...<i.d.+i(a;), (2.25) 

rr I X ( dyVd^...dada+Ax) \'^ , d'^yVd^...dsds+Ax) , c roORA 

Ud....dXx) = [ j + + (2.26) 

These determine a new Hamiltonian system with s + 1 deletions: 

Ud, ...da^, = Ad, ...da^.Al ,„,^^^ + Sd^^, = + U,i, ... (x), (2.27) 

J J def ( dWd,...ds+M V d''Wd,...ds+,{x) , ^ . . 

f/.i (a:) - j + ^..,., (2.28) 

</'di...(i,+in(a:;) =^^(ii...d,+i0di...d.n(a;) (n G . . . , 4+i}), (2.29) 

'Hdi...d,+i(pdi...d,+inix) ^ Sn(f)di...ds+in{x) {u G Z+\{c/i, . . . , 4+i}). 



(2.30) 



After deleting all the T> = {di, ■ ■ ■ , di} energy levels, the resulting Hamiltonian system 
= 'Hdi...dc: Av = Adi...di-, etc has the following form: 

Uv = AvA^ + Ed, = + Uv{x), (2.31) 
^^^'^ " V — — / — ^2?(3;) = log(/>di...d,(a;), (2.32) 

(^Vn{x)^^ Av(^dx-dt-^n{x) {u & Z+\D) , (2.33) 
^T5<^Pn(a;) = £:n075„(x) (nGZ+\P). (2.34) 
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Now that "Hd has the lowest energy level /i: 

^ = mm{n\neZ+\V}, (2.35) 

with the groundstate wavefunction 0^(x) 

(pf,{x) = (I)t:>^{x) = 0di...d,^(x). (2.36) 

As usual the Hamiltonian system can be expressed simply in terms of the groundstate wave- 
function (j)fj.{x), which we will denote by new symbols 7/, A, etc: 

'H = nv = A^A + = p^ + U{x), (2.37) 
Uix) ^ U^,{x) (^^)' + + ^.' ^(x) ^ W^,(x) '='log0,(x), (2.39) 

mn{x)=£nK{x). ^n{x) = (t)Vn{x) {u E Z+\V) . (2.40) 

As shown by Krein-Adler [18], the results can be expressed succinctly: 

Mx) = -T777- — 7 , X ('^ e Z+\P), (2.41) 

W[0d,,0d2, . . . ,0dj(a;) 

^(:^^) = f/d,...A(a;) = U{x) - 2^(logW [(^d,, . . . , (l)dM^)) > 0), (2.42) 



in which the Wronskian determinant is defined by 

'd^~'fkixy^ 

l<j,k<n 

For n = 0, we set W [■]{x) = 1. In deriving the determinant formulas f l2.4ip and fl2.42p use 
is made of the properties of the Wronskian 



W[/....,/„](.)^^^det(^) . (2.43) 



W[gh,gf2, . . . , gfn]{x) = (7(x)"W[/i, h,..., /„](x), (2.44) 

W[W[/i,/2,...,/n,^7],W[/i,/2,...,/n,/i]](x) 

= W[/i, /2, . . . , fn]{x) W[/i, /2, ...,fn,g, h]{x) (n > 0). (2.45) 

Let us note that Udi...de{x) and (pdi,...den{x) are symmetric with respect to di,...,di, and 
thus Ti = 1-Ldi ...di is independent of the order of {dj}. 

Let us state Adler's theorem again; // the set of deleted energy levels di, . . . ,di satisfy 
the condition (12. Sp . the Hamiltonian 1-i = 'Hdi...di = p'^ + U{x) with (12.420 is well-defined 
and hermitian, and its complete set of eigenf unctions (I2.40p are given by fl2.4ip . Crum's 
theorem corresponds to the choice {di, . . . , de} = {0, 1, — 1}, and the resulting lowest 
energy level is fj, = i and there is no deleted energy levels above the new groundstate. 



2.2 Polynomial eigenf unctions 

In this subsection we consider the typical case of shape invariant systems in which the 
eigenfunctions consist of the orthogonal polynomials {Pn}'- 

= M^)PrM^)), M^) = e'^("\ (2.46) 

in which t]{x) is called the sinusoidal coordinate. As shown in detail in the examples in 
Appendix A, 77(2;) = x for the harmonic oscillator (the Hermite polynomials), r]{x) = x^ for 
the radial oscillator (the Laguerre polynomials) and ri{x) = cos2x for the Darboux-Poschl- 
Teller potential (the Jacobi polynomials). The groundstate wavefunction <po{x) provides the 
orthogonality weight function 

Ipoixf Pn{r]{x))Pm{ri{x))dx = hnSnm, U, 171 E Z+. (2.47) 



XI 



In this case, the modification of Crum's theorem produces the eigenfunctions {0„(x)} 
which again consist of polynomials in 'r]{x). By using fl2.44l) and 

/,(x)'=V,(r/(a;)), W[A, Z^, /„](x) = (^) ^"^""Vi/i, /s, ...,/„] (77(0;)), (2.48) 

we obtain a simple expression of the eigenfunctions 

M^) - M^) (^) w[p.„p,„...,p.j(,(x)) • ^'-''^ 

This simply means that the resulting eigenfunctions are again polynomials in r]{x): 

Mx)=iP{x)Vniv{x)), (2.50) 

^(-'" w|^.t'S(.M) ' ^~M"W|P..P......P..P„](,). (2.51) 

satisfying the orthogonality relation 



X2 
XI 



'>P{xf 'Pn{ri{x))Vm{'n{x))dx = KSnm, Tl, 171 E Z+\V. (2.52) 



Let us emphasise that n is not the degree in i] and by construction, £ members are missing: 
Vdi = Vd2 = ■ ■ ■ = = 0. Therefore these polynomials cannot be called exceptional 
orthogonal polynomials |El [121 IT5] . 
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3 'Discrete' Quantum Mechanics 



Let us begin with a few general remarks on the one-dimensional discrete QM with pure 
imaginary shifts. See |25] for the general introduction to the discrete quantum mechanics 
with pure imaginary shifts and [34] for the Crum's theorem in the discrete QM. In the discrete 
QM, the dynamical variables are, as in the ordinary QM, the coordinate x, which takes value 
in an infinite or a semi-infinite or a finite range of the real axis and the canonical momentum p, 
which is realised as a differential operator p = —id^. Since the momentum operator appears 
in exponentiated forms e''^'^^, 7 G M, in a Hamiltonian, it causes finite pure imaginary shifts 
in the wavefunction e'^'^^ip{x) = ip{x =F ^7). This requires the wavefunction as well as other 
functions appearing in the Hamiltonian to be analytic in x within a certain complex domain 
including the physical region of the coordinate. Let us introduce the ^-operation on an 
analytic function, * : / /*. If /(x) = ^ e C, then f*{x) = E a*x", in which a* 

n n 

is the complex conjugation of an- Obviously = f{x) and f{x)* = f*{x*). If a function 

satisfies /* = /, we call it a 'real' function, for it takes real values on the real line. 

The starting point is again a generic one dimensional discrete quantum mechanical Hamil- 
tonian with discrete semi- infinite energy levels only (12. ip . Again we assume that the ground- 
state energy is chosen to be zero Sq = 0, so that the Hamiltonian is positive semi-definite. 
The generic factorised Hamiltonian reads [221 125] 

n = A^A = y^V{x) e''P^yV*{x) + e'^P^/V{x) - V{x) - V*{x), (3.1) 

A=U{e^P./[^-e-^P^/V{^), A^ = -t{^/V{^e^P-y^V^e-iP). (3.2) 

Since the ^-operation for Af, A'^ f and l-if satisfies 

{AfYix) = Arix), (A^fYix) = A^rix), wnx) = nnx), (3.3) 

they map a 'real' function to a 'real' function 

r = f {Afr = Af, {A^fr = A^f, {nfr = nf. (3.4) 

By specifying the function V{x), various explicit examples are obtained [221 [251 E]. A few 
exactly solvable examples are given in Appendix. The corresponding Schrodinger equation 

Hi/j^x) = Si/j^x) is a difference equation 

\/V{x)V*{x — ij) ip{x — ij) + ^/V*{x)V{x + i'j) ip{x + ^7) 
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- {V{x) + V*{x))'iIj{x) = £iJj{x), (3.5) 

instead of differential in the ordinary QM. Although this equation looks rather compli- 
cated, the equation for the polynomial eigenfunctions (13.471) has a familiar form of difference 
equations. Again the groundstate wavefunction (j)Q{x) is determined as a zero mode of A, 
A4>o{x) = (^ 'H(j)o{x) = 0), namely, 



,Jv*{x - i^)Mx - i^) - ^Jv{x + t^)Mx + ^i) = 0. (3.6) 

This dictates how the 'phase' of the potential function V is related to that of the groundstate 
wavefunction 0o- Here we also assume that the groundstate wavefunction (f)o{x) has no node 
and chosen to be real and positive for real x. 

Due to the lack of generic theorems in the theory of difference equations, let us assume 
that all the energy levels are non-degenerate and that all the eigenfunctions are square- 
integrable and orthogonal with each other and form a complete basis of the Hilbert space: 

'H(f)n{x) = £n4>n{x), TL G Z+, (3.7) 

/ (t)n{x)*(j)m{x)dx = hn^nmj < /l„ < CX), 71,171 (3.8) 
J Xl 

In most explicit examples these statements can be verified straightforwardly. For simplicity 
we choose all the eigenfunctions to be real on the real axis 0* = 0n, which is made possible 
by dH. 

3.1 Modification of Crum's theorem 

The formulation of the modified Crum's theorem in the discrete quantum mechanics goes 
almost parallel to that in the ordinary quantum mechanics. Again the presentation is purely 
algebraic. Let us note that various quantities in intermediate steps might have singularities 
and Hamiltonians might not be hermitian. We choose a set of distinct non-negative integers 
V =^ {di,d2, ■ ■ ■ ,de} C Z+, satisfying the condition (12.81) as before. First let us note that 
the Hamiltonian T-L can be rewritten by incorporating the level di as: 

n = AlAd, + Sd„ Ad,(t>dr = 0, (3.9) 
^,/=%(ei^^/i^-e-i^v/T^), 4/=^-^(v'T^e^^-y^V^e-^^), (3.10) 

(3.11) 



VM ^V{x)V*{x-^^) "^'fj^^^^ ^ 
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These define a new Hamiltonian system 

nd,^^Ad,A\^+£d,, (3.12) 

y'di4>din{x) = Sn(pdin{x), (pdi n{x) = Adr4>n{x) , U & Z+\{di} . (3.13) 

Note that the energy level di is now deleted, (j)^^ di (x) = 0, from the set of 'eigenf unctions' 
{<Pdin{x)} of the new Hamiltonian T-Ldi- 

Suppose we have determined T-Ldi...ds and 0di ...d^ ra(a^) with s deletions. They have the 
following properties 

..... = Ad, ...d Ai...d,+S d^, (3.14) 
Ad, ..... = ^ (e^^^olj^ - e-^^WVd,...dM), 

^i...d. = -i{VVd,...dAx) e^^ - ^Jvi___,^{x) e-^^), (3.15) 
|Vn.)^*(x-n)^^^ (. = 1), 

(Pdi...dsn{x) = Ad,...dAdi...d,-,n{x), 4>d,...dsn{x) = (t)^ ...dsui^) ^ (3-17) 

T^d,...ds4>di...dsn{x) — £n4>di ...dsn{x), (3.18) 

where n e Z+\{(ii, . . . , dg}. We have also 

(t>di...d,-in{x) ^ j-^^^^(f>di...dsn{x) {neZ+\{di,...,ds}). (3.19) 

Next wc will define a new Hamiltonian system with one more deletion of the level dg+i- We 
can show the following: 

T-ldi...ds = '^d^...dsds+1'^d.i...dsds+i + ^ds+i, Ad,...dsds+i'Pdi...dsds+i{x) ^0, (3.20) 

= -^(VK.,...d...,,(x) ei- - ^Vl,,,,^,^Jx) e-i^), (3.21) 

K.. ......... (x) JVd,...d.{x - zi)V-,:...,.(x - zi) '^^-••^-^-^-^'^ "'"^^ (3.22) 

V <?>.i...d.(i.+i(,a:;j 

These determine a new Hamiltonian system with s + 1 deletions: 

'Hdi...ds+i = >^di ...d.+ivAj^i + ^d.+i) (3.23) 
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(/}d,...d,+,n{x)=Ad,...d,+,(l)d,...d,n{x), (l)d,...d.,+,n{x) = <P*d, ...d^+^ni^) , (3.24) 

'Hdi ...d^+i4'di ...ds+in{x) = £n4>di...ds+in{x)-, (3.25) 

where n G Z+\{(ii, . . . , ds+i}- 

After deleting all the V = {di, ■ ■ ■ , di} energy levels, the resulting Hamiltonian system 
T-Lt) = 'Hdi...di, Av = •Adi...di, etc has the following form: 

Hv = AvAl + Sd,, (3.26) 

Av=U(^e^P^/l^)-e-^^^/VU^), ^J, = -^(v/T^ei^ - ^l^e-i^) , (3.27) 

Vv{x) = ^Vd, ...d,.A^ - il)Vl ...,^_^(x - .i) , (3.28) 

4>Vn{x)=Av(t)d^...dt_^n{x), (pv n{x) = (p^ ni^) , {u E Z+\V) , (3.29) 

T-iv(pVn{x) = £n4>Vn{x) {u G '^+\D). (3.30) 

Now that T/x) has the lowest energy level /i: 

/i =^ min{n I n G Z+yP}, (3.31) 

with the groundstate wavefunction <j)fj,{x) 

(j)^{x) =^ (pvfiix) = 0di ■■■d,/x(x). (3.32) 

Then the Hamiltonian system can be expressed simply in terms of the groundstate wave- 
function (pfxix), which we will denote by new symbols "H, A, etc: 

H = nv = A^A + Sf,, A^^{x) = 0, (3.33) 



A = Av, = I {ei^sjv*{x) - e-^^y^) , 

V{x) ^ y^^(x) = JVv{x - zi)F^(x - il) ^^^"^ ~ , (3.35) 

'H4>n{x) = £n4>n{x), (^„(x) = (pVn{x) {u G Z+\P). (3.36) 

The discrete counterpart of the determinant formulas ( I2.4ip - (]2.42p requires a deformation 
of the Wronskian, the Casorati determinant, which has a good limiting property: 

W,[/i, . . . , Q{x) det(/fc(x + 2^^7)) , (3.37) 

\ / l<j,k<n 
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lim 7-^"("-')W,[/i,/2, . . ■Jn]{x) = W[/i,/2, . . . ,/n](x), 
7^0 



(3.38) 



(for n = 0, we set [■](^) = !•)• satisfies 

W,[/i,...,/„]*(a;)=W,[/*,...,/:](a;), 

n 

W,[gh,gf2, ...,gfn] = l[ g{x + ^^^7) " W,[/i, /a, . . . , /„](x), 

[W^[/l, /2, . . . , /n, ^7], W^[/i, /2, . . . , /„, h] ] (X) 

= W^[/i, /2, . . . , /„](x) W^[/i, /2, . . . , /n, ^7, h]{x) (n > 0). 



(3.39) 
(3.40) 

(3.41) 



By using the Casorati determinant we obtain (£ > 0) 



0„(x) = (t)di...dtn{x) 



-i-i 



7) 



[0di,...,0dJ(a;-«^)' 



(3.42) 



y(x) = = ^T/(x-2|7)T/*(x-z^7) 

[0rf,,...,0rfj(x + if) [(/)d,,...,(/)rf,,0^](a:-z7) 



, . . . , 0dj {x - a) , . . . ,(|)d„(p^l] (x) 



(3.43) 



We also have {£ > 0) 



,,,(x + ^^7) 



1 



j=0 



(3.44) 



Therefore Vrfj...rf^(a;) and (/>di,...d^n(a:) are symmetric with respect to di, . . . ,di, and Tidi^.d,, 
is independent of the order of {dj}. 

Let us state the discrete QM analogue of Adler's theorem; // the set of deleted energy 
levels V = {di, . . . ,di} satisfy the condition f l2.8p . the modified Hamiltonian is given by 
"H = 'Hdi...dt = A^A + S^ with the potential function given by (I3.43P and its eigenf unctions 
are given by fl3.42p . The discrete QM version of Crum's theorem [31] corresponds to the 
choice {(ii, . . . , di} = {0, 1, 1} and the new groundstate is at the level fi = i and 

there is no vacant energy level above that. Due to the lack of generic theorems in the 
theory of difference equations, the hermiticity of the resulting Hamiltonian "H and the non- 
singularity of the eigenfunctions (f>n{x) cannot be proved categorically for the discrete QM, 
even when the condition (12. 8p is satisfied by the deleted levels. See Appendix A of [25] for a 
detailed discussion of the self-adjointness of the Hamiltonians in discrete QM. It should be 
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stressed that in most practical cases, in particular, in the cases of polynomial eigenfunctions, 
the hermiticity of the Hamiltonian 7/ and non-singularity of the eigenfunctions are 
satisfied. 

3.2 Polynomial eigenfunctions 

In this subsection we consider the typical case of shape invariant systems in which the 
eigenfunctions consist of the orthogonal polynomials {Pn}'- 

Mx) = Mx)Pn{v{x)), AMx) = 0, (3.45) 

in which //(x) is called the sinusoidal coordinate, see (IB. 91) . The groundstate wavefunction 
(f)o{x) provides the orthogonality weight function 



(poixY Pn{ri{x))Pm{,'n{x))dx = hnSnm, U, TTl & Z+. (3.46) 

The difference equation for {Pn} looks much simpler than the Schrodinger equation (13.51) : 

V{X) {Pnivix - 27)) - Pnivix))) + V* (x) {Pn{v{x + l^)) " Pn{v{x))) = EnPuHx)). (3.47) 



For the explicit forms of V{x), see for example ( IB. 101) . these are the equations that determine 



the hypergeometric orthogonal polynomials, e.g. the Meixner-Pollaczek (MP), the continu- 
ous Hahn (cH), the Wilson (W) and the Askey- Wilson (AW) polynomials. In fact, the above 
form of the difference equation fl3.47p is independent of the fact that is a polynomial or 
not. It is obtained simply by the similarity transformation of the Hamiltonian (13. ip in terms 
of the groundstate wavefunction (po^x): 



n = Mx)~^ ° ^ ° Mx) = V{x) + V{x)* e-^P - V{x) - V{xy. (3.48) 

In the case of polynomial eigenfunctions, the modification of Crum's theorem produces 
the eigenfunctions {0„(x)} which again consist of polynomials in t]{x). By using the property 
fl3:i0|l we have 

e 

W^ [</>!,... , (Pf]{x) = J](/)o(x + t'-^i) ■ W^[A, . . . , Pi]{x), (3.49) 
i=i 

i+i 

W^[0i, ...,(/),, <j)n\{x) = l[<j)o{x + I'-^i) ■ W^[A, ....h, Pn]{x), (3.50) 
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where Pn{x) == Pn{f]{x)). Corresponding to the formula fl2.48p in the ordinary QM, we have 



fj{x) = fjivi^)), fjiv)- polynomial in r], 

W^[/i, /2, . . . , fn]{x) = (pn{x) X (polynomial in r]{x)), 



{3.51] 



in which ipn{x) is defined in flB.39P (and Lp{x) is defined in flB.20p ). 

These simply mean that the resulting eigenfunctions {(/)n{x)} (I3.42p are again polynomials 
in 77(x): 



0„(x) = i/j^x) Vn{ri{x)), 



(3.52) 



def 



. WVd^,„d^{x + i 



^+1-,- . ^i+i{x) 0o(a^ + %7) 
2 ^^^^_q)Q^r,{x-iT)) 



\ l[^{x-^|^)^V{x + ^'-^^)■ 
\ k=0 



0o(x - i^j) 

Qivix-tW 



X 



ll'fi(x + il~f)Jv{x-i^fi',') 



k=0 



Q(r/(x + a)) ' 



(3.53) 



in which Vn{v{x)) and 2(77(2;)) are certain polynomials in 77(2;) defined by 



W^[Frf,, . . . , Pd„Pn] = ^e+iix) X Vnivix)), W^[Frf,, . . . , = ipeix) x Qivix)). (3.54) 

The polynomials {Vn} form a complete basis of the Hilbert space and satisfy the orthogo- 
nality relations 

rX2 



ip{x) Vni'nix))Vmiv{x))dx = KSnm, n, m G Z+\P. 



(3.55) 



XI 



Let us emphasise that n is the level of the original eigenfunction and not the degree in r/. 
The degree of Vnijj) depends on the set and it can be calculated explicitly from f l3.54p . 
By construction t members are missing: Vd^ = Vd2 = ■ ■ ■ = Vd^ = 0. Therefore these 
polynomials cannot be called exceptional orthogonal polynomials [9l [121 115] . 



4 Summary and Comments 

Theory of exactly solvable discrete QM is less developed than that of the ordinary QM. 
Up to date, the known exactly solvable discrete quantum systems are all shape invariant 
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[Ml [25] and in one to one correspondence with the known (g)-hypergeometric orthogonal 
polynomials [26]-|29]. Small progress was made in this direction |31] by introducing several 
new sinusoidal coordinates for the construction of new types of exactly solvable Hamiltonians. 
Roughly speaking, this approach attempts to create the discrete analogues of various Morse 
type potentials and the soliton potentials. In this paper we pursue another direction; to 
construct infinitely many exactly solvable quantum systems by deforming the known exactly 
solvable one. In the ordinary QM, the modification of Crum's theorem [I] due to Krein-Adler 
[18] allows to produce an essentially iso-spectral Hamiltonian by deleting a finite number of 
energy levels from the original system. The set of deleted level must satisfy certain condition 
(12. 8p . but there are infinitely many possible deletions leading to infinitely many exactly 
solvable systems starting from a known one. The discrete analogue of Adler's modification 
is presented in this paper in parallel with the original version, since the algebraic structure 
is common. We also comment on the practical cases when the eigenfunctions consist of 
orthogonal polynomials. The eigenfunctions of the resulting system also consist of orthogonal 
polynomials. But certain members of these polynomials are missing due to the deletion. 

Very special and simple examples, in which all the excited states from the first to the i-th 
are deleted (see Fig. 2), are presented explicitly in Appendix. As will be commented shortly, 
these examples were instrumental for the discovery of the infinitely many shape invariant 
systems and the corresponding infinitely many exceptional orthogonal polynomials [121 [T3] . 
In the ordinary QM, the corresponding prepotential has a very simple form flA.lSp : 

w,{x- A) = W{x- A + M) + log \ 

to be compared with the prepotential for the Hamiltonian of the ^-th exceptional orthogonal 
polynomial (14) and (28) of [12]: 

w,{x; A) = Wix; X + i6) + log ^^(^(^)^ ^ + ^) . (4.1) 

In the discrete QM, the corresponding formula is ( IB. 351) 

Ve{x; X) = K —— .. V{x; X + iS), 

^e{r]{x-%j);X) 

to be compared with the corresponding formula for the Hamiltonian of the £-th exceptional 
orthogonal polynomial (30) of [T3] : 

Ve{x; A = — . f f \ \ < x\ — TTl ^TVTT^^' ^ + M ■ 4.2 

^i{ri{xy, X + 5) ^e{r]{x - i^); X) 
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The addition (multiplication) of the deforming polynomial with the shifted parameters 
^e{ri{x)] \ + 6) would achieve the shape invariance. Since the harmonic oscillator has no 
shiftable parameter, we have ^£{rj{x); A) = ^e{v{^)'j A + 5). This also 'explains' non-existence 
of exceptional Hermite polynomials. In contrast to the Hermite polynomial, the continuous 
Hahn polynomial has four real parameters. We can construct the corresponding exceptional 
continuous Hahn polynomials with three real parameters, which will be reported elsewhere. 

The actual function forms of the deforming polynomial in Appendix are not the same 
as those for the exceptional orthogonal polynomials. For the ordinary QM examples, see 
( 1X23]) vs. (13) and (27) in [E] and for the discrete QM examples, see (iRiTj) vs. (64) and 
(78) in |13]. But they share some interesting features. 

Before closing this section, let us remark that the present modification of Crum's theorem 
is applicable to the Hamiltonian systems of various species of the infinite family of exceptional 
orthogonal polynomials p^[T3l[T5] . as well as to those of the classical orthogonal polynomials 
including the Wilson and the Askey- Wilson polynomials. 
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Appendix 

In Appendix we present very special and simple examples of an application of Adler's theo- 
rem, in which the eigenstates 0i, 02, • • • , 0^ are deleted. In other words, V = {di, d2, . . . , de}= 
{1,2,...,^}, that is, the modified groundstate level is the same as that of the original theory 
/i = 0. The situation is illustrated in Fig. 2, which should be compared with Fig. 1, depicting 
the generic case discussed in sections 12.11 and 13.11 
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Figure 1 : Generic case 



Figure 2: Special case 



The black circles denote the energy levels, whereas the white circles denote deleted energy 
levels. We write H = 7/12. 0n = 0i2...^n, A = A12...1, V = Vv2...i etc. as Ui, At, 
etc. This Hamiltonian = AIAi is non-singular for even i but may be singular for odd i. 
Since algebraic formulas such as the Wronskians and Casoratians are valid for even and odd 
i, we present various formulas without restricting to the even i. The original systems are 
shape invariant but the . . . , 0£)-deleted systems He are not. The rightmost vertical line 
in Fig. 2 corresponds to the Hamiltonian system "H^ = AtA\, which is shape invariant and it 
is obtained from l-Ln by one more step of Crum's method. This study helped us to find the 
new shape invariant systems and exceptional orthogonal polynomials 



A The ordinary QM 

Here we apply Adler's theorem to the harmonic oscillator, the radial oscillator and the 
Darboux-Poschl- Teller potential, whose eigenfunctions are described by the classical orthog- 
onal polynomials. That is, the Hermite, Laguerre and Jacobi polynomials, to be abbreviated 
as H, L and J, respectively. These original systems are shape invariant, meaning a very spe- 
cial form of parameter dependence, f lA.Sp . flA.Op . Here we display the parameter dependence 
explicitly by A, which represents the set of the parameters. 
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A.l The original systems 



Here we summarise various properties of the original Hamiltonian systems to be compared 
with the specially modified systems to be presented in IA.2I Let us start with the Hamilto- 
nians, Schrodinger equations and eigenf unctions {xi < x < X2): 



def ,,,^+ . / x n def (i>V(x; A) _ rf>V(x; A) 



?/(A) ^=^U(A)U(A), ^(A)^='- 



dx 



dx 



n{X)Mx; A) = SnWMx; A) (n = o, i, 2, . . .), 

A) = 0o(x; A)P„(?7(x); A), M^', A) = e^(^^^). 



(A.l) 
(A.2) 
(A.3) 



Here iji^x) is the sinusoidal coordinate, W(x; A) is the prepotential and £^n(A) is the n-th 
energy eigenvalue: 



ri{x) 



def 





Xi 


= —00, 


X2 


= 00, 


: H 




Xi 


= 0, 


X2 


= 00, 


: L 


COS 2x, 


Xi 


= 0, 


X2 


TV 

~ 2' 


: J 



dcf 



A) 



def 



1^2 
2 

■\x'^ + glogx 



H 

L , £n{X) 

J 



def 



none 

9, ^ > 
(51, /i), g,h>0 

2n 
An 

An{n + g + h) 



H 

L , (A.4) 
J 



H 

L 
J 



(A.5) 



g log sinx + h log cos x 

The eigenfunction consists of an orthogonal polynomial Pn{v] A), a polynomial of degree n 
in r), {Pniv'j A) = for n < 0): 



i^.(r];A)'='c„(A)P~(r;;A), 

Hniv) ■■ H 



L , Cn{X) 
J 



def j (-1)" 
n! 

(n+g+fc) 



2"n! 



in which (a)„ is the Pochhammer symbol. Shape invariance means 



Ai\)A{X)^ = A{\ + SyAi\ + d)+Si{X), 6={ 1 



none 



or equivalently, 



5)\2 d^W{x;X + d) 



H 
L 
J 



H 
L 
J 



(A.6) 
(A.7) 



(A.8) 



dx"^ 



+ ^i(A). 



(A.9) 



/ dWjx; A) \2 _ rf2>v(^. A) _ A 
\ dx J dx"^ V 

The action of ^(A) and A{\y on the eigenfunction is: 

^(A)0,(x;A) = /„(A)(/*,_i(x;A + (5), ^(A) Vn-i (x; A + 5) = 6„_i(A)(^,(x; A). (A.IO) 
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Here the coefficients /n(A) and 6„_i(A) are the factors of Sn{X): 



2n 



-2{n + g + h) 



H 



L, ^n-iW = 1^^271 -LJ' ^n(A) = /„(A)6„_i(A). 

(A.ll) 

The forward and backward shift operators, J^(A) and B{X), are defined by: 

J-(A) Mx; A + d)-' o ^(A) o A) = J'^'':^., ^ , (A.12) 

(po{x; \ + o) ax 

BiX) = 0o(x; A)-^ o A{Xy o ,^o(x; X + 6) 

and their action on the polynomial is: 

J^{X)PMx); A) = UX)Pn-i{v{xy, X + 6), (A.14) 
BiX)P,_,{r]{x)- X + d) = b^_^{X)PMxy, A). (A.15) 

Note that -7^(A) and B{X) can also be expressed in terms of r] only [TB]. The orthogonality 
reads 



X2 



o(x; A) Pnivix); X)Pm{,r]{x)] X)dx = /i„(A)5„„, (A. 16) 



{2"n!y^ : H 

^r(n + ^? + i) :L. (A.17) 

r(n+g+|)r(n+fc+|) ^ , 
2n\{2n+g+h)r(n+g+h) 

A. 2 The (01, . . . , 0£)-deleted systems 

The prepotential of the modified system is obtained from (I2.5ip up to an additive constant: 

we{x; A) = log0,,o(x) = W(x; X + £6) - hgUvi^); A). (A.18) 

It is a polynomial {^i{ri{x); A)) deformation of the shape invariant one W(x; A + £6). Note 
that the normalization of does not affect the Hamiltonian. The explicit forms of the 
deforming polynomial ^^{r]; A) will be given in ( ]A.23p . For even £, the polynomial ^£{r]{x); A) 



has no zero in the range of x and the modified Hamiltonian system is hermitian, which reads: 

dcf dwejx; A) ^ d dwejx; A) 
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7{,(A) = A(A)U,(A), 

Hi{\)<Pt,n{x- A) = f„(A)0,,„(x; A) (n = 0, £ + 1, £ + 2, . . .) 



(A.20) 
(A.21) 



This system is not shape invariant. As mentioned in section HJ the above form of the 
deformed prepotential f lA.lSp is closely related to that of the exceptional Laguerre and Jacobi 
polynomials. 

A degree I polynomial in r], ^^(z/; A) is defined by 



W[Fi, . . . , P,]{r^- A) = n k\ Cfc(A) ■ A), 
and the explicit forms are: 



(A.22) 



fe=i 



Urn A) 



r (-9-^-1) 



{g+h+l)i ^1 



-ri) 



iv) 



H 
L 
J 



(A.23) 



The eigenfunctions are 

(pe,o{x] A) = e 



def Wf{x;X) 



bo{x; A + id) 



e 

W[Pi, . . . , Pe, PnKv, A) = n ^! c,(A) ■ {-lYPeAv, A) Pi,o{v; A) = l) • 



(A.24) 
(A.25) 



k=l 



Note that Pi^niv'^ A) is a polynomial of degree nmrj and Po,n{v'^ A) = Pn(?7; A) and Pi^niVi A) = 
for 1 < 72 < We set Pe^niv'i A) = for n < 0. For even £, the eigenpolynomial Pg^nivi^)'! A) 
(n > £ + 1) has n — £ zeros in the range of x. The operators ^^(A) and ^^(A)^ connect the 
modified system 'Hf(A) = Ai{\y Ae{X) to the shape invariant system 'H^(A) = ^£(A)^^(A)''" 
with the parameters X + {£ + 1)6, which is denoted by the rightmost vertical line in Fig. 2. 
The n-th level (n > £ + 1) of the modified system Tie is iso-spectral with the n — £ — 1-th 
level of the new shape invariant system T-L'f 



Ae{X)(j)e,n{x] X) = /£,„(A)0„_£_i (x; X + {£+ 1)5) , 
A(A)Vn-^-i [x] X+{£ + 1)5) = 6^,n-i(A)0£,„(x; A). 



(A.26) 
(A.27) 



Here fi^nW and be^n-iW are defined by 



An(A) = /n(A)x A 



be,n-iiX) = bn-i{X)xA-\ 



A 



[-2)\n-£)e 
I 

:-2)-^(n + (7 + /i + l), 



H 
L 
J 



(A.28) 
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and they factorise £^n(A), £n{^) = f£,nWbe,n-i{^) ■ The forward and backward shift opera- 
tors J^i{X) and Be{X), which act on the polynomial eigenfunctions, are defined by: 



J-,(A) = (Po{x; X + {e+ 1)6) ' o MX) o (^,,0(2;; A) 
_ (j)o{x; X + eS) 1 d 

~ (/.o(x;A + (£ + l)5)e£(r/(a;);A) ' 

Bi{X) 0,,o(2:; A)-i o ^,(A)t o ,^o(x; A + + 1)6) 



-1 



(A.29) 



{x;X+{i+l)6) 



[x- X + M) 



iein{x)\X) 



dx^iivi^), ^) 



X I — + 4(>V(x;A + M) + >V(x;A + (£+l)(5)) 



(A.30) 



Their action on the polynomials is {n > £ + 1): 



TeiX)Pe,Mx); X) = /,,„(A)P„_,_i (vix); X+{i+ 1)6) , 
BiiX)Pn-i-i{vix); A + (£ + 1)6) = h,n-i{X)PMx); A). 



(A.31) 
(A.32) 



Note that J-'^(A) and -^^(A) can be expressed in terms of i] [16]. For n > I + 1, the above 
relation (IA.32P provides a simple formula of the modified eigenpolynomial Pe^niv'i ^) i'^ terms 
of (,i{T]; X) and the original eigenpolynomial PniVi ^)- 

btn-liX)fn-i{X + i6)Pi^n{v;X) 

= £n-iiX + i6)Uv; X)Pn-t{m a + M) + Ac2{7])dMm A) d^Pn-t{m A + M), (A.33) 



in which the coefficient 02(77) is given by 



C2{v) 



def 



H 
L 
J 



(A.34) 



The orthogonality relation for even i is: 



>; A)^ Pi^nivix); X)Pl^rn{'n{x)] X)dx = hi^n{X)Snm, 



XI 



he,n{X) = {n - i)e hn{X) x <( 1 

4~^{n + g + h + l)i 



H 

L , (n = 0, n > £ + 1). 
J 



(A.35) 
(A.36) 



A few historical remarks are in order. Dubov et al [T7] derived in 1992 an exactly solvable 
Hamiltonian system of a deformed harmonic oscillator, which corresponds to the i = 2 case 
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of this Appendix. Their paper, written about two years before Adler's, rehed on rather 
heuristic arguments. Recently Quesne [11] derived exactly solvable and non-shape invariant 
systems of deformed radial oscillator and deformed DPT potential, both are called type 
III. These results again correspond to the ^ = 2 cases of the radial oscillator and the DPT 
potential of this Appendix. 

B The discrete QM 

Here we apply Adler's theorem to the shape invariant, therefore solvable, systems whose 
eigenf unctions are described by the orthogonal polynomials; the Meixner-Pollaczek (we set 
the parameter = ^), continuous Hahn, Wilson and Askey- Wilson polynomials [29J, to be 
abbreviated as MP, cH, W and AW, respectively. See [22] and [25] for the discrete QM 
treatment of these polynomials. 

B.l The original systems 

Here we summarise various properties of the original Hamiltonian systems to be compared 
with the specially modified systems to be presented in IB.2I Let us start with the Hamilto- 
nians, Schrodinger equations and eigenfunctions {xi < x < X2): 

A{\) ^=^%(eiV^*(^; A) - e-"rPy/V{x- A) ), 

^(A)t ^=^f -^{^/V{x■, A) - y/V*{x; A) e'^^'), (B.l) 

?^(A) = ^(A)U(A), (B.2) 

7^(A)(/.„(x;A) =^„(A)</)„(a;;A) (n = 0, 1, 2, . . .), (B.3) 

0„(x; A) = 0o(x; X)Pn{vix); A). (B.4) 

The set of parameters A are 

MP: X = a, a>0, (B.5) 

cH: X={ai,a2), Rea^ > (i = l,2), (B.6) 
W: A =^ (ai, 02, as, 04), Reoj > {i = 1, . . . ,4), 

{a*, Cg, ttg, 04} = {ai, 02, (23, 04} (as a set), (B.7) 
AW: g-^ =^ (ai, 02, as, 04), |ai| < 1, (i = l,...,4), 0<q<l, 

{a^, ttg, Og, 04} = {ai, 02, as, 04} (as a set), (B.8) 
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where g('*^i'-^2,.- ) {q'^^, g^^, . . .). The the sinusoidal coordinate r]{x) is, 



?7(x) 



dcf 



X, 
„2 



Xi = — OO, X2 = oo, 7=1 
X", Xi = 0, X2 = OO, 7=1 

cosx, Xi = 0, X2 = vr, 7 = logq' 



The potential function V{x; A) and energy eigenvalue ^„(A) are 



V{x;X) = { 



' a + ix 
(fli + 2a;) (a2 + ix) 



IX] 



def 



(22x(2ix + l))" ni=i(«j + 
^ ((1 - e2-)(l - ge^-))-' n'=i(l - %e-) 

2n 



MP, cH 

W 

AW 



MP 

cH 

W 
AW 



(B.9) 



(B.IO) 



def 



6i = ai + a2 + a* + 

def 



(B.ll) 



n{n + 6i - 1), 

n(n + 6i - 1), 6i = ai + a2 + a3 + 04 

(g"" - 1)(1 - M""^), ^4 =^ 01020304 

The eigenfunction is described by the orthogonal polynomial Pnijl'i ^ polynomial of degree 
n in rj: 



MP 
cH 
W 
AW 



P„(r,;A)'=l'c„(A)Pr"'^(r/;A), 



4^f J Pn(?7; 01,02,01,02) 
Wn{r); 01,02,03,04) 
p„(?7;oi,02,03,04|g) 



^n(?7;A) 



MP 
cH 
W 
AW 



^2" 

77.! 



Cn(A) 



def 



^(n + &i-l)„ 
(-l)"(n + 6i-l)„ 

2"(64g"-^g)n 



MP 
cH 
W 
AW 



(B.12) 



, (B.13) 



in which (a; g)„ is the g-Pochhammer symbol. We set P„(?7; A) = for n < 0. The shape 
invariance relations involve one more parameter k: 



A{X)A{Xy = kA{\ + 5)U(A + d)+ £i{X), 

MP 



def 



'1 i i 1^ 
^2' 2' 2' 2> 



cH , K. 

W,AW 



def I 1 : MP,cH,W 



^-1 :AW 



(B.14) 
(B.15) 



or equivalently. 



V{x - q; X)V*{x - q; A) = V{x; A + d)V*{x - 27; A + 5), 
V{x + q; A) + V*{x- q; A) = k{V{x; X + d) + V*{x; X + d))- Si{X). 

The groundstate wavcfunction (f)o{x]X) is determined by 

^V*{x-q;X) 00 - ^i; A) = ^V{x + q;X) 0o(x + A), 



(B.16) 
(B.17) 

(B.18) 
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and its explicit forms are: 



(x; A) 



def 



A/r(ai + ix)T{a2 + ix)T{al — ix)T{a2 — ix) 
{V{2ix)V{-2ix))-^ n5=i r(aj- + ix)T{aj - ix) 



We introduce an auxiliary function (p{x) with the properties: 



MP 
cH 

W 

AW 



(B.19) 



ip{x] 



def 



1 

2x 

2 sin X 



MP, cH 

W 

AW 



(f)o{x; A + d) = (p{x)^V{x + A) 0o(a: + A), 



(p{x) 

The sinusoidal coordinate r]{x) has the following properties: 



ri{x - i^) - r){x + i^) = -i(p{x) x 



k 



: MP, cH,W 



r]{x - i^) + r]{x + i^) 



r]{x-i^)r]{x + i^) 



2ri{x) 

2ri{x) - 
2?7(x) cosh ^ 
r]{xY + 

?7(a;)2 + sinh^ ^ 



sinh^ :AW 



MP, cH 

W 

AW 

MP, cH 

W 

AW 



(B.20) 

(B.21) 
(B.22) 

(B.23) 
(B.24) 

(B.25) 



These mean that for a polynomial P{r]) in rj, i(p(x) ^[P(r)(x — i^)) — P(r)(x + i^))) is 
another polynomial in r]{x). The action of ^(A) and A{Xy on the eigenfunctions is 

A{X)Mx;X)^fn{X)(j>n-i{x;X + S), ^(A)Vn-i(a:; A + 6) = 6„_i(A)<^„(x; A). (B.26) 

The factors of the energy eigenvalue, /n(A) and bn-i{X), Sn{X) = /n(A)6n-i(A), are given 

by 

(B.27) 



/ , s def J n + 6i — 1 
JniAj — 



—n{n + 6i — 1) 
gf (g-ri - 1)(1 - 64^"-!) 



MP 
cH 
W 
AW 



, bn-i{X)''^'{ -1 

n 

q 2 



MP, cH 
W 
AW 



The forward and backward shift operators J-'{X) and ^(A) are defined by 
J^(A) = 0o(^; A + 6)-^ o ^(A) o (j)o{x; A) = iip{x)-\e^P - e'^P), 



(B.28) 
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B{X) = (Po{x; X)-^ o ^(A)^ o (j)o{x; X + d) = -i{V{x; X)e^P - V*{x; X)e-2P)^{x), (B.29) 



and their action on the polynomials is 



J-(A)F„(r/(a;); A) = /„,(A)P„_i(?7(x); X + 5), 
B{X)Pn-iH^)-, X + 5) = hn-i{X)PnHx)- A). 



The orthogonality relation is 



X2 



KiX) 



A)^ Pn{r]{x)] A)P„(77(x); X)dx = hn{X)5nm, 

' 27rr(n + 2a)(n!22«)~^ 
dcf j 27rn?j=i r(n + a, + a*) ■ (n!(2n + 61 - \)Y{n + bi - 1))"^ 
27m! {n + bi-l)nU 

i<j<j<4 -'^('^ + fli + c^j) ■ r(2r?, + 61) ""^ 
[ 2n{b,q-''; qUb.q'''; q)oo{?^'; q)^' ni<.<,<4(«^%9"; q)^' 



MP 
cH 
W 
AW 



(B.30) 
(B.31) 



(B.32) 



(B.33) 



B.2 The (01, 



-deleted systems 



The potential function, the Hamiltonian and the Schrodinger equation of the modified system 
are: 



dei ^{x - i^-f--f)ip{x - i^-f) ^e{r]{x + qy, X) 



^e{r]{x - q); X) 
MX) = 2(e^yy;(x; A) - e-i*'^A^(^[A) ) , 



^,(A)t 'i^* -^(v/\4(x;A)e^^' - A) e~^P) , 

7/,(A) = A(A)U,(A), 

?^^(A)0^,n(a;; A) = ^„(A)0^,„(x; A) (n = 0, £ + 1, ^ + 2, 



(B.34) 
(B.35) 

(B.36) 
(B.37) 
(B.38) 



The explicit forms of the deforming polynomial ^^(r^; A) will be given in f lB.41|) . For even i, 
C,e{r]{x)] A) has no zero in the rectangular domain in the complex x plane, xi < Rex < X2, 
|Imx| < ||7|. Note that the normalization of does not affect Tie. This system is not shape 
invariant. The second line of the expression for Ve{x] A), flB.35p . is obtained from the first 
line (1B.34P by repeated use of the formula (1B.22P of the auxiliary function 99. As mentioned 
in section HI this form of the deformed potential function ( ]B.35p is closely related to that of 
the exceptional Wilson and Askey- Wilson polynomials. 
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It is convenient to introduce an auxiliary function tpeix): 

£-2 

Mx) = Vixp m^ix - 2l^Mx + (B.39) 

k=l 

where [x] denotes the greatest integer not exceeding x. Note that [|]+2 Ylk=il^T^] ~ |^(^— !)• 
The deforming polynomial ^eirj; A) is defined by 

W^[A, . . . , mx; A) A) = Pnivix); A)) 

Ul=ikl :MP,cH,W 



def 

k=l 

As in the ordinary QM cases ( ]A.23p . it is expressed in terms of the polynomial Pe of the 
original system with shifted parameters: 

P,{r^--y-{i-l)S) f jiiy' :MP,cH,W 

^) = c,i-y-ii-i)S) " 1 iUU -F)- : AW • ^^-^^^ 

Note that F„(r^; A*) = Pn{v; A) and c„(A*) = c„(A) for the MP, W and AW cases. The 
eigenfunctions are 

(pefi{x;X} = (pe^nix; X) = (pe,o{x;X)Pe,n[V{x);X), (B.42) 



\^^eivix - ?f ); X)^i{r]{x + ); A) 
W^[Pi, . . . , P^, Pn]{x; A) = J] Cfc(A) ■ ^e+i{x){-lYPiAv{xy, A) (B.43) 

A;=l 

' nLi^! : MP,cH,W 



nLin;=isinh^ :aw 

For even £, Pi^n{v{x)', A) {n > i + 1) has n — £ zeros in the range of x. Note that Pe^niv'i ^) 
is a polynomial of degree tj, in and Fo,n(^; A) = Pniv'i ^) P£,n{v'i A) = for 1 < n < £. 
We set Pe^niv'i A) = for n < 0. The operators ^^(A) and ^^(A)''" connect the modified 
system 'H^(A) = ^^(A)"''^£(A) to the shape invariant system '^^(A) = ^^(A)^^(A)"I' with the 
parameters X+ {£ + 1)6, which is denoted by the rightmost vertical line in Fig. 2. The n-th 
level [n > £ + 1) of the modified system Tie is iso-spectral with the n — £ — 1-th level of the 
new shape invariant system l-i'f 

At{X)(t)t,n{x; X) = fe,niX)^n-i-i {x; X+i£+ 1)6) , (B.44) 
AiiX)^(j)n-i-i{x;X+i£ + l)6) =h,n-i{X)(Pi,n{x;X). (B.45) 
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Here, fe,nW and 6^_„_i(A) are the factors of the energy eigenvalue, £n{^) = fe,nWbe,n-iW 
and are defined by 

' 2' 



An(A) = /„(A)x A J ih + n 



i-lYin - iUh + n)e 



MP 

W • 

AW 



The forward and backward shift operators J^e{X) and Bi{X) which act on the polynomial 
eigenf unctions, are defined by: 

TeiX) = 00 (x; A + (£ + 1)6)-' o Ae{X) o </),,o(x; A) = ±-^^^^{e'^P - e'^^), (B.47) 

BiiX) = <Pe,o{x; Xy' o A,(A)t o 0o(x; A + (£ + 1)6) 

= ^-iY^-W-^)^-i) (v{x; X + i6)Uvix + ^i); A)eiP 

- V*{x; X + £6)Uvix - ^i); A)e-^p) (B.48) 
and their action on the polynomials is 

T,{X)Pe,n{v{x);X) = fe,n{X)Pn-e-i{v{x);X + {e+l)6), (B.49) 
BiiX)Pn-e-i{vix);X+ie + l)6) = 6,,„_i(A)P,,„(?7(x); A). (B.50) 

For n > i + 1, the above formula (IB. 501) provides a simple formula of the modified eigen- 
polynomial Pg^niv'^ A) in terms of CiiV'j A) and the original eigenpolynomial Pn{v] X): 

{-lYK\'^'-'\n-l{X)P,,n{m A) 

= -i{y{x- X + (i6)Uri{x + il)- X)ip{x - z|)P„_,_i(r7(x - q)- A + + 1)5) 

- V*{x- X + i6)Uv{x - tl)] X)^{x + «|)F„_,_i(?7(x + 2|); A + (£ + 1)(5)) . (B.51) 

The orthogonality relation for even (. is: 

I (j)i>fl{x] X)^ Pe^nivix); X)Pi^^{r]{x); X)dx = hi^niX)5nm, (B.52) 

(n - i)e2^ : MP 

he^^{X) = K{X) X \ {n- i)eibi + n)e : cH,W , (n = 0,n >£ + !). (B.53) 

AW 



g ^"(g« ^■,q)i{hq"';q)i 
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